Solar Physics 
DOI: 10.1007/' 



Magnetic Connectivity between Active Regions 
10987, 10988, and 10989 by Means of Nonlinear 
Force-Free Field Extrapolation 



Tilaye Tadesse^ ^, T. Wiegelmann^, and 
B. Inhester^, and A. Pevtsov^ 



© Springer •••• 



Abstract Extrapolation codes for modelling the magnetic field in the corona in 
cartesian geometry do not take the curvature of the Sun's surface into account 
and can only be applied to relatively small areas, e.g., a single active region. We 
apply a method for nonlinear force-free coronal magnetic field modelling of pho- 
tosphcric vector magnctograms in spherical geometry which allows us to study 
the connectivity between multi-active regions. We use vector magnctograph data 
from the Synoptic Optical Long-term Investigations of the Sun survey (SO- 
LIS)/Vector Spectromagnetograph(VSM) to model the coronal magnetic field, 
where we study three neighbouring magnetically connected active regions (ARs: 
10987, 10988, 10989) observed on 28, 29, and 30 March 2008, respectively We 
compare the magnetic field topologies and the magnetic energy densities and 
study the connectivities between the active rcgions(ARs). We have studied the 
time evolution of magnetic field over the period of three days and found no major 
changes in topologies as there was no major eruption event. From this study we 
have concluded that active regions are much more connected magnetically than 
the electric current. 
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1. Introduction 

In order to model and understand the physical mechanisms underlying the var- 
ious activity phenomena that can be observed in the solar atmosphere, like, for 
instance, the onset of fiarcs and coronal mass ejections, the stability of active 
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region, and to monitor the magnetic helicity and free magnetic energy, the mag- 
netic field vector throughout the atmosphere must be known. However, routine 
measurements of the solar magnetic field are mainly carried out in the pho- 
tosphere. The magnetic field in the photosphere is measured using the Zeeman 
effect of magnetically-sensitive solar spectral lines. The problem of measuring the 
coronal field and its embedded electrical currents thus leads us to use numerical 
modelling to infer the field strength in the higher layers of the solar atmosphere 
from the measured photospheric field. Except in eruptions, the magnetic field 
in the solar corona evolves slowly as it responds to changes in the surface field, 
implying that the electromagnetic Lorentz forces in this low-/3 environment are 
relatively weak and that any electrical currents that exist must be essentially 
parallel or antiparallel to the magnetic field wherever the field is not negligible. 

Due to the low value of the plasma /? (the ratio of gas pressure to magnetic 
pressure), the solar corona is magnetically dominated ( [Gary, 2001[ ). To describe 
the equilibrium structure of the static coronal magnetic field when non-magnetic 
forces are negligible, the force-free assumption is appropriate: 



(V X B) X B = 



(1) 



V-B = 



(2) 



B^^g on photosphere 



(3) 



where B is the magnetic field and Bp|-,g is measured vector field on the photo- 
sphere. Equation ((TJ states that the Lorentz force vanishes (as a consequence of 
J II B, where J is the electric current density) and Equation (0) describes the 
absence of magnetic monopoles. 

The extrapolation methods based on this assumption are termed nonlinear 



force- free field extrapolation ( Sakurai, 1981 Amari et ai, 1997 Amari, Boulmezaoud, and Mikic, 1999H 



Amari, Boulmezaoud, and Aly, 2006 Wu et al, 1990 Cuperman, Demoulin, and Semel, 1991 1 



Demoulin, Cuperman, and Semel, 1992 Inhester and Wicgelmann, 2006 Mikic and McClymont, 1994H 



Roumeliotis, 1996 Yan and Sakurai, 2000 Valori, Kliem, and Keppcns, 2005 



Wheatland, 2004^, Wheatland and Rcgnier, 2009[ [Wheatland and Lcka, 2010 



Wiegelmann, 2004 1 



Amari and Aly, 2010] ) .| 



For a more complete review of existing methods for computing nonlinear force- 
free coronal magnetic fields, we refer to the review papers bv lAmari et flLI(jl997p . 
[Schrijver ef"ar]((MM)) . [Metcalf et aZipM)) . and |WiegeliiIaSI| ([innHl). | Wiegelmann and Neukirch| 
(|2006p has developed a code for the self-consistent computation of the coronal 
magnetic fields and the coronal plasma that uses non-force-free MHD equilibria. 

The magnetic field is not force-free in the photosphere, but becomes force- 
free roughly 400 km above the photosphere ( [Metcalf et ai, 1995[ ). Furthermore, 
measurement errors, in particular for the transverse field components (i.e., per- 
pendicular to the line of sight of the observer), would destroy the compati- 
bility of a magnetogram with the condition of being force-free. One way to 
ease these problems is to preprocess the magnetogram data as suggested by 
[Wiegelmann, Inhester, and Sakur"ai|(j2006( ). The preprocessing modifies the bound-| 
ary values of B within the error margins of the measurement in such a way that 
the moduli of force-free integral constraints of Molodcnskii ( [Molodensky, 1974| ) 
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are minimized. The resulting boundary values are expected to be more suitable 
for an extrapolation into a force-free field than the original values. 

In the present work, wc use a larger computational domain which accom- 
modates most of the connectivity within the coronal region. We also take the 
uncertainties of measurements in vector magnetograms into account as suggested 
in IDeRosa et al.\ (|2009p . We apply a preprocessing procedure to SOLIS data 
in spherical geometry ( Tadesse, Wiegelmann, and Inhester, 2009D by taking ac- 
count of the curvature of the Sun's surface. For our observations, performed on 
28, 29, and 30 March 2008, respectively, the large field of view contains three 
active regions (ARs: 10987, 10988, 10989). 

The full inversion of SOLIS/VSM magnetograms yields the magnetic filling 
factor for each pixel, and it also corrects for magneto-optical effects in the 
spectral line formation. The full inversion is performed in the framework of 
Milne-Eddington model (ME)( [Unno, 1956] ) only for pixels whose polarization 
is above a selected threshold. Pixels with polarization below threshold are left 
undetermined. These data gaps represent a major difficulty for existing magnetic 
field extrapolation schemes. Due to the large area of missing data in the example 
treated here in this work, the reconstructed field model obtained must be treated 
with some caution. It is very likely that the field strength in the area of missing 
data was small because the inversion procedure, which calculates the surface 
field from the Stokes line spectra, abandons the calculation if the signal is below 
a certain threshold. The magnetic field in the corona, however, is dominated by 
the strongest flux elements on the surface, even if they occupy only a small 
portion of the surface. We are therefore confident that these dominant flux 
elements are accounted for in the surface magnetogram, so that the resulting 
field model is fairly realistic. At any rate, it is the field close to the real, which 
can be constructed from the available sparse data. Therefore, we use a procedure 
which allows us to incorporate measurement error and treat regions with lacking 
observational data as in lTadesse et a/.l(|201ip . The technique has been tested in 
cartesian geometry in [Wiegelmann and Inhester] (|2010p for synthetic boundary 
data. 



2. Optimization Principle in Spherical Geometry 



Wheatland, Sturrock, and Roumeliotis (j2000p have proposed the variational prin-| 



ciple to be solved iteratively which minimizes Lorentz forces ([T]) and the diver- 
gence of magnetic field ^ throughout the volume of interest, V. Later on the 
procedure has been improved by |Wiegelmann| (|2004p for cartesian geometry in 
a such way that it can only uses the bottom boundary on the photosphere as an 
input. Here we use optimization approach for the functional (Ci^) in spherical ge- 
ometry ( [Wiegelmann, 2007] [Tadesse, Wiegelmann, and Inhester, 2009P and iter- 
ate B to minimize C^j. The modification concerns the input bottom boundary 
field Bpi-jg which the model field B is not forced to match exactly but wc allow 
deviations of the order of the observational errors. The modified variational 
problem is ( [Wiegelmann and Inhester, 2010[ [Tadesse et ai, 201ip : 

B = argmin(£(^) 
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C^ = C^ + C^ + i^^photo (4) 

£f= / u;f(r, 6*, 1 (V X B) X BjV^ sin 6'(ird6'd(/) 

Jv 

Jv 

Wo= / (B - Bobs) •W(0,0). (B-Bobs)r2 sin 
^ s 

where and C^-^ measure how well the force-free Equations ([T]) and divergence- 
free ^ conditions are fulfilled, respectively. a;£(r, 0, 0) and ^, are weight- 
ing functions for force-free term and divergence-free term, respectively and iden- 
tical for this study. The third integral, i^^ihoiO' ^ surface integral over the pho- 
tosphere which serves to relax the field on the photosphere towards a force-free 
solution without too much deviation from the original surface field data, B^bg- In 
this integral, W(0, </>) = diag('u;].^jjg^b '^trans' ^trans) diagonal matrix which 
gives different weights for observed surface field components depending on its 
relative accuracy in measurement. In this sense, lacking data is considered most 
inaccurate and is taken account of by setting W{9, (f) to zero in all elements of 
the matrix. 

We use a spherical grid r, 9, (j) with ri^, ng, grid points in the direction 
of radius, latitude, and longitude, respectively. In the code, we normalize the 
magnetic field with the average radial magnetic field on the photosphere and 
the length scale with a solar radius for numerical reason. Figure [T] shows a map 
of the radial component of the field as color-coded and the transverse magnetic 
field depicted as white arrows. For this particular dataset, about 86% of the data 
pixels are undetermined. The method works as follows: 

• We compute an initial source surface potential field in the computational 
domain from B^bg tbe normal component of the surface field at the pho- 
tosphere at r = IRq. The computation is performed by assuming that a cur- 
rentless (J = or V x B = 0) approximation holds between the photosphere 
and some spherical surface Ss (source surface where the magnetic field 
vector is assumed radial) . We compute the solution of this boundary- value 
problem in a standard form of spherical harmonics expansion. 

• We minimize C^^ (Equations |4|) iteratively. The model magnetic field B at 
the surface is gradually driven towards the observed field B^bg while the 
field in the volume V relaxes to force-free. If the observed field, B^bg: 
is inconsistent, the difference B — B^bg remains finite depending in the 
control parameter v. At data gaps in B^bg: the respective field value is 
automatically ignored. 

• The iteration stops when becomes stationary as t^L^j < 10^"', 
is the decrease of Ccj during an iterative steps. 

• A convergence to = yields a perfect force-free and divergence-free state 
and exact agreement of the boundary values B with observations B^bg in 
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-1200 -960 -720 -480 -240 240 480 720 960 1200 
Figure 1. Surface contour plot of radial magnetic field component and vector field plot 
of transverse field with white arrows. The color coding shows Br on the photosphere. The 
vertical and horizontal axes show latitude, 8{in degree) and longitude, (/'(in degree) on the 
photosphere. In the area coloured in olive, field values are lacking. The region inside the black 
box corresponds to the physical domain where the weighting function is unity and the outside 
region is the buffer zone where it declines to zero. The blue boxes indicate the domains of the 
three active regions. 



regions where the elements of W are greater than zero. For inconsistent 
boundary data the force-free and solenoidal conditions can stiU be fulfilled, 
but the surface term ^pY^Q^Q will remain finite. This results in some de- 
viation of the bottom boundary data from the observations, especially in 
regions where ^^.g^j^jg^j and wtrans small. The parameter v is tuned so 
that these deviations do not exceed the local estimated measurement error. 



3. Results 

In this work, we apply our extrapolation scheme to Milne-Eddington inverted 
vector magnetograph data from the Synoptic Optical Long-term Investigations 
of the Sun survey (SOLIS). As a first step for our work we remove non-magnetic 
forces from the observed surface magnetic field using our spherical preprocessing 
procedure. The code takes B^^g as improved boundary condition. 

SOLIS/VSM provides full-disk vector- magnetograms, but for some individ- 
ual pixels the inversion from line profiles to field values may not have been 
successful inverted and field data there will be missing for these pixels (see 
Figure [1]). The different errors for the radial and transverse components of B^j^g 
are taken account by different values for Wj-^^j^j and u^trans • ^^^^ work wc used 
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Figure 2. Wedge-shaped computational box of volume V with the inner physical domain V' 
and a buffer zone. O is the center of the Sun. 



^radial ~ ^'^'^trans ^^'^ surface preproccsscd fields as the radial component 
of B^i^g is measured with higher accuracy. 

We compute the 3D magnetic field above the observed surface region inside 
wedge-shaped computational box of volume V, which includes an inner phys- 
ical domain V' and a buffer zonc(thc region outside the physical domain), as 
shown in Figure [5J The physical domain V' is a wedge-shaped volume, with 
two latitudinal boundaries at ^j^^jj^ = —26° and ^rnax = 16°, two longitudinal 
boundaries at (/Jmin ~ 129° and ipmax = 226°, and two radial boundaries at 
the photosphere (r = IRq) and r = 1.75i?0. Note that the longitude cj) is 
measured from the center meridian of the back side of the disk. We define 
V' to be the inner region of V (including the photospheric boundary) with 
w-f = u!^ ~ 1 everywhere including its six inner boundaries SV'. We use a 
position-dependent weighting function to introduce a buffer boundary of nd = 10 
grid points towards the side and top boundaries of the computational box, V . 
The weighting functions, and uj^ are chosen to be unity within the inner 
physical domain V' and decline to with a cosine profile in the buffer bound- 
ary region ( Wiegelmann, 2004 Tadesse, Wiegelmann, and Inhester, 2009[ ). The 



framed region in Figure [T] corresponds to the lower boundary of the physical 
domain V' with a resolution of 114 x 251 pixels in the photosphere. 

The middle panel of Figure [3] shows magnetic field line plots for three con- 
secutive dates of observation. The top and bottom panels of Figure |3] show the 
position of the three active regions on the solar disk both for SOLIS full-disk 
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(a) 28 March 2008 15:45UT (b) 29 March 2008 15:48UT (c) 30 March 2008 15:47UT 






(d) 28 March 2008 15:45UT (c) 29 March 2008 15:48UT (f) 30 March 2008 15:47UT 



(g) 28 March 2008 16:00UT (h) 29 March 2008 15:48UT 



(i) 30 March 2008 15:48UT 



Figure 3. Top row: SOLIS/VSM magnetograms of respective dates. Middle row: Magnetic 
field lines reconstructed from magnetograms on the top panel. Bottom row: EIT image of the 
Sun at I95A on indicated dates. 



magnetogram Q and SOHO/EItI image of the Sun observed at 195A on the 
indicated dates and times. Figure 0] shows some selected magnetic field lines 
from reconstruction from the SOLIS magnetograms, zoomed in from the middle 
panels of Figure |3l In each column of Figure |4] the field lines are plotted from the 
same foot points to compare the change in topology of the magnetic field over the 
period of the three days of observation. In order to compare the fields at the three 
consecutive days quantitatively, we computed the vector correlations between the 
three field configurations. The vector correlation (Cvec) ( jSchrijver et ai, 2006[ ) 
metric generalizes the standard correlation coefficient for scalar functions and is 



^http:/ /solis. nso.edu/solis_data. html 

^http: / /sohowww. nascom.nasa.gov/data/archive 
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(c) 30 March 2008 15:47UT 



Figure 4. Some magnetic field fines plots reconstructed from SOLIS magnctograms using 
nonlinear force-free modelling. The color coding shows Br on the photosphere. 
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Table 1. The magnetic energy associated with extrapolated NLFF 
field configurations for the three particular dates. 



Date 


EnlffClO^^erg) 


Epot(1032erg) 


£;free(1032erg) 


28 March 2008 


57.34 


53.89 


3.45 


29 March 2008 


57.48 


54.07 


3.41 


30 March 2008 


57.37 


53.93 


3.44 



given by 

where and arc 3D vectors at grid point i. If the vector fields are identical, 
then Cvec = 1; if _L , then CVec = 0. The correlations (Cycc) of the 3D 
magnetic field vectors of 28 and 30 March with respect to the field on 29 March 
are 0.96 and 0.93 respectively. From these values we can see that there has been 
no major change in the magnetic field configuration during this period. 

We also compute the values of the free magnetic energy estimated from the 
excess energy of the extrapolated field beyond the potential field satisfying 
the same B^^^t, ■ f boundary condition. Similar estimates have been made by 
[Regnier and Priest] (|2007p and [Thalmann, Wiegelmann, and Raouafi] (|2008p for 
single active regions observed at other times. From the corresponding potential 
and force- free magnetic field, BpQ^; and B, respectively, we can estimate an 
upper limit to the free magnetic energy associated with coronal currents 

i?frec - i?nlfl[ - ^pot = -tJy, (^nlff - ^pot) r^^^^OdrdOd^. (6) 



The computed energy values are listed in Table [TJ The free energy on all three 
days is about 3.5 x lO'^^erg. The magnetic energy associated with the potential 
field configuration is about 54 x lO'^^erg. Hence E^^Q exceeds i?pot by only 6%. 
Figure [5] shows Iso-surface plots of magnetic energy density in the volume above 
the active regions. There are strong energy concentrations above each active 
region. There were no major changes in the magnetic energy density over the 
observation period and there was no major eruptive phenomenon during those 
three days in the region observed. 



In our previous work(Tadesse et al., 20111, we have studied the connectivity 



between two neighbouring active regions. In this work with an even larger field 
of view, the three ARs share a decent amount of magnetic flux compared to their 
internal flux from one polarity to the other (see Figure |3]). In terms of the electric 
current they are much more isolated. In order to quantify these connectivities, 
we have calculated the magnetic flux and the electric currents shared between 
active regions. For the magnetic flux, e.g., we use 

= ^ |B, • f|i?|sin(0,)A0,A(/), (7) 
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(a) 28 March 2008 15:45UT 




(b) 29 March 2008 15:48UT 




(c) 30 March 2008 15:47UT 



Figure 5. Iso-surfaccs (ISs) of the absolute NLFF magnetic energy density (7.5 X lO^^erg) 
for the three consecutive dates computed within the entire computational domain. 
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Table 2. The percentage of the total magnetic flux shared between the three ARs. $11 , ^22 
and ^33 denote magnetic flux of AR 10989{lcft), AR 10988{middle) and AR 10987(right) 
of Figure^ respectively. 

28th 29th 30th 

= 1 2 3 a = l 2 3 a = 1 2 3 

/3 = 1 56.37 5^59 OOO 56.50 5^48 OOO 56.50 5^48 0.00 

2 13.66 81.12 1.43 13.66 81.22 1.43 13.66 81.22 2.22 

3 0.00 0.48 71.47 0.00 0.48 71.80 0.00 0.48 71.80 
Elsewhere 29.97 12.82 27.10 29.84 12.82 26.77 29.84 12.82 25.98 



Table 3. The percentage of the total electric current shared between the three ARs. 
^11, h2, and 733 denote electric current of AR 10989(left), AR 10988(middle) and AR 
10987(right) of Figure [T] respectively. 

28th 29th 3Qth 

1^/3 = 1 2 3 a = 1 2 3 a = 1 2 3 

/3 = 1 82.47 019 OOO 86.36 019 OOO 94.16 019 0.00 

2 0.65 85.25 1.42 0.65 85.25 1.42 0.65 85.25 3.55 

3 0.00 0.38 82.27 0.00 0.38 82.27 0.00 0.38 82.27 
Elsewhere 16.88 14.18 16.31 12.99 14.18 16.31 5.19 14.18 14.18 



where the summation is over all pixels of ARq from which the field line ends in 
AR^ or i G ARqII conjugate footpoint(i) € ARp. The indices a and /3 denote the 
active regions and the index number 1 corresponds to AR 10989, 2 to AR 10988 
and 3 to AR 10987 of Figure [T] For the electric current we replace the magnetic 
field, Bi • f, by the vertical current density J,; • f in Equation ([?])• Whenever 
the end point of a field line falls outside (blue rectangles in Figure [T]) of the 
three ARs, we categorize it as ending elsewhere. Both Table [2] and [3] show the 
percentage of the total magnetic fiux and electric current shared between the 
three ARs. So, for example first column of Tabled] shows that 56.37% of positive 
polarity of ARi is connected to negative polarity of ARi; line 2 shows that 
13.66% of positive/negative polarity of ARi is connected to positive/negative 
polarity of AR2, and line 3 shows that there are no field lines (0%) connecting 
positive/negative polarity of ARi with positive/negative polarity of AR3. The 
same technique applies for Table |3]too. The three active regions are magnetically 
connected but much less by electric currents. 



4. Conclusions 



We have investigated the coronal magnetic field associated with three ARs 10987, 
10987, 10989, on 28, 29 and 30 March 2008 by analysing SOLIS/VSM data. We 
have used an optimization method for the reconstruction of nonlinear force-free 
coronal magnetic fields in spherical geometry by restricting the code to limited 
parts of the Sun ( Wiegelmann, 2007| [Tadesse, Wiegelmann, and Inhester, 2009[ 



Tadesse et ai, 20111. The code was modified so that it allows us to deal with 



lacking data and regions with poor signal-to-noise ratio in a systematic manner( Wiegelmann and Inhester, 201 
[Tadesse et al, 2011] ). 



SOLA: solar_paper.tex; 19 January 2013; 12:37; p. 11 



T. Tadesse et al. 



We have studied the time evolution of magnetic field over the period of three 
days and found no major changes in topologies as there was no major eruption 
event. The magnetic energies calculated in the large wedge-shaped computa- 
tional box above the three active regions were not far apart in value. This is 
the first study which contains three well separated ARs in our model. This was 
made possible by the use of spherical coordinates and it allows us to analyse 
linkage between the ARs. The active regions share a decent amount of magnetic 
flux compared to their internal flux from one polarity to the other. In terms of 
the electric current they are much more isolated. 
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